\(\int \frac {d+e x^2+f x^4}{x^4 (a+b x^2+c x^4)} \, dx\) [59]

Optimal result

Integrand size = 30, antiderivative size = 267 \[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=-\frac {d}{3 a x^3}+\frac {b d-a e}{a^2 x}+\frac {\sqrt {c} \left (b d-a e+\frac {b^2 d-a b e-2 a (c d-a f)}{\sqrt {b^2-4 a c}}\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a^2 \sqrt {b-\sqrt {b^2-4 a c}}}-\frac {\sqrt {c} \left (b^2 d-b \left (\sqrt {b^2-4 a c} d+a e\right )-a \left (2 c d-\sqrt {b^2-4 a c} e-2 a f\right )\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a^2 \sqrt {b^2-4 a c} \sqrt {b+\sqrt {b^2-4 a c}}} \]

[Out]

-1/3*d/a/x^3+(-a*e+b*d)/a^2/x+1/2*arctan(x*2^(1/2)*c^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2))*c^(1/2)*(b*d-a*e+(b^2
*d-a*b*e-2*a*(-a*f+c*d))/(-4*a*c+b^2)^(1/2))/a^2*2^(1/2)/(b-(-4*a*c+b^2)^(1/2))^(1/2)-1/2*arctan(x*2^(1/2)*c^(
1/2)/(b+(-4*a*c+b^2)^(1/2))^(1/2))*c^(1/2)*(b^2*d-b*(a*e+d*(-4*a*c+b^2)^(1/2))-a*(2*c*d-2*a*f-e*(-4*a*c+b^2)^(
1/2)))/a^2*2^(1/2)/(-4*a*c+b^2)^(1/2)/(b+(-4*a*c+b^2)^(1/2))^(1/2)

Rubi [A] (verified)

Time = 0.68 (sec) , antiderivative size = 267, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 3, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.100, Rules used = {1678, 1180, 211} \[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=\frac {\sqrt {c} \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b-\sqrt {b^2-4 a c}}}\right ) \left (\frac {-a b e-2 a (c d-a f)+b^2 d}{\sqrt {b^2-4 a c}}-a e+b d\right )}{\sqrt {2} a^2 \sqrt {b-\sqrt {b^2-4 a c}}}-\frac {\sqrt {c} \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {\sqrt {b^2-4 a c}+b}}\right ) \left (-a \left (-e \sqrt {b^2-4 a c}-2 a f+2 c d\right )-b \left (d \sqrt {b^2-4 a c}+a e\right )+b^2 d\right )}{\sqrt {2} a^2 \sqrt {b^2-4 a c} \sqrt {\sqrt {b^2-4 a c}+b}}+\frac {b d-a e}{a^2 x}-\frac {d}{3 a x^3} \]

[In]

Int[(d + e*x^2 + f*x^4)/(x^4*(a + b*x^2 + c*x^4)),x]

[Out]

-1/3*d/(a*x^3) + (b*d - a*e)/(a^2*x) + (Sqrt[c]*(b*d - a*e + (b^2*d - a*b*e - 2*a*(c*d - a*f))/Sqrt[b^2 - 4*a*
c])*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*a^2*Sqrt[b - Sqrt[b^2 - 4*a*c]]) - (Sqrt
[c]*(b^2*d - b*(Sqrt[b^2 - 4*a*c]*d + a*e) - a*(2*c*d - Sqrt[b^2 - 4*a*c]*e - 2*a*f))*ArcTan[(Sqrt[2]*Sqrt[c]*
x)/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[2]*a^2*Sqrt[b^2 - 4*a*c]*Sqrt[b + Sqrt[b^2 - 4*a*c]])

Rule 211

Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[a/b, 2]/a)*ArcTan[x/Rt[a/b, 2]], x] /; FreeQ[{a, b}, x]
&& PosQ[a/b]

Rule 1180

Int[((d_) + (e_.)*(x_)^2)/((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4), x_Symbol] :> With[{q = Rt[b^2 - 4*a*c, 2]}, Di
st[e/2 + (2*c*d - b*e)/(2*q), Int[1/(b/2 - q/2 + c*x^2), x], x] + Dist[e/2 - (2*c*d - b*e)/(2*q), Int[1/(b/2 +
q/2 + c*x^2), x], x]] /; FreeQ[{a, b, c, d, e}, x] && NeQ[b^2 - 4*a*c, 0] && NeQ[c*d^2 - a*e^2, 0] && PosQ[b^
2 - 4*a*c]

Rule 1678

Int[(Pq_)*((d_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^2 + (c_.)*(x_)^4)^(p_.), x_Symbol] :> Int[ExpandIntegrand[(d*x
)^m*Pq*(a + b*x^2 + c*x^4)^p, x], x] /; FreeQ[{a, b, c, d, m}, x] && PolyQ[Pq, x^2] && IGtQ[p, -2]

Rubi steps \begin{align*} \text {integral}& = \int \left (\frac {d}{a x^4}+\frac {-b d+a e}{a^2 x^2}+\frac {b^2 d-a b e-a (c d-a f)+c (b d-a e) x^2}{a^2 \left (a+b x^2+c x^4\right )}\right ) \, dx \\ & = -\frac {d}{3 a x^3}+\frac {b d-a e}{a^2 x}+\frac {\int \frac {b^2 d-a b e-a (c d-a f)+c (b d-a e) x^2}{a+b x^2+c x^4} \, dx}{a^2} \\ & = -\frac {d}{3 a x^3}+\frac {b d-a e}{a^2 x}+\frac {\left (c \left (b d-a e-\frac {b^2 d-a b e-2 a (c d-a f)}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\frac {b}{2}+\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx}{2 a^2}+\frac {\left (c \left (b d-a e+\frac {b^2 d-a b e-2 a (c d-a f)}{\sqrt {b^2-4 a c}}\right )\right ) \int \frac {1}{\frac {b}{2}-\frac {1}{2} \sqrt {b^2-4 a c}+c x^2} \, dx}{2 a^2} \\ & = -\frac {d}{3 a x^3}+\frac {b d-a e}{a^2 x}+\frac {\sqrt {c} \left (b d-a e+\frac {b^2 d-a b e-2 a (c d-a f)}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a^2 \sqrt {b-\sqrt {b^2-4 a c}}}+\frac {\sqrt {c} \left (b d-a e-\frac {b^2 d-a b e-2 a (c d-a f)}{\sqrt {b^2-4 a c}}\right ) \tan ^{-1}\left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\sqrt {2} a^2 \sqrt {b+\sqrt {b^2-4 a c}}} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.21 (sec) , antiderivative size = 284, normalized size of antiderivative = 1.06 \[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=\frac {-\frac {2 a d}{x^3}+\frac {6 b d-6 a e}{x}+\frac {3 \sqrt {2} \sqrt {c} \left (b^2 d+b \left (\sqrt {b^2-4 a c} d-a e\right )+a \left (-2 c d-\sqrt {b^2-4 a c} e+2 a f\right )\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b-\sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {b-\sqrt {b^2-4 a c}}}+\frac {3 \sqrt {2} \sqrt {c} \left (-b^2 d+b \left (\sqrt {b^2-4 a c} d+a e\right )-a \left (-2 c d+\sqrt {b^2-4 a c} e+2 a f\right )\right ) \arctan \left (\frac {\sqrt {2} \sqrt {c} x}{\sqrt {b+\sqrt {b^2-4 a c}}}\right )}{\sqrt {b^2-4 a c} \sqrt {b+\sqrt {b^2-4 a c}}}}{6 a^2} \]

[In]

Integrate[(d + e*x^2 + f*x^4)/(x^4*(a + b*x^2 + c*x^4)),x]

[Out]

((-2*a*d)/x^3 + (6*b*d - 6*a*e)/x + (3*Sqrt[2]*Sqrt[c]*(b^2*d + b*(Sqrt[b^2 - 4*a*c]*d - a*e) + a*(-2*c*d - Sq
rt[b^2 - 4*a*c]*e + 2*a*f))*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b - Sqrt[b^2 - 4*a*c]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[b
- Sqrt[b^2 - 4*a*c]]) + (3*Sqrt[2]*Sqrt[c]*(-(b^2*d) + b*(Sqrt[b^2 - 4*a*c]*d + a*e) - a*(-2*c*d + Sqrt[b^2 -
4*a*c]*e + 2*a*f))*ArcTan[(Sqrt[2]*Sqrt[c]*x)/Sqrt[b + Sqrt[b^2 - 4*a*c]]])/(Sqrt[b^2 - 4*a*c]*Sqrt[b + Sqrt[
b^2 - 4*a*c]]))/(6*a^2)

Maple [A] (verified)

Time = 0.12 (sec) , antiderivative size = 244, normalized size of antiderivative = 0.91

method	result	size
default	\(-\frac {d}{3 a \,x^{3}}-\frac {a e -b d}{a^{2} x}+\frac {4 c \left (\frac {\left (-a e \sqrt {-4 a c +b^{2}}+b d \sqrt {-4 a c +b^{2}}-2 f \,a^{2}+a b e +2 a c d -b^{2} d \right ) \sqrt {2}\, \arctan \left (\frac {c x \sqrt {2}}{\sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) c}}\right )}{8 \sqrt {-4 a c +b^{2}}\, \sqrt {\left (b +\sqrt {-4 a c +b^{2}}\right ) c}}-\frac {\left (-a e \sqrt {-4 a c +b^{2}}+b d \sqrt {-4 a c +b^{2}}+2 f \,a^{2}-a b e -2 a c d +b^{2} d \right ) \sqrt {2}\, \operatorname {arctanh}\left (\frac {c x \sqrt {2}}{\sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) c}}\right )}{8 \sqrt {-4 a c +b^{2}}\, \sqrt {\left (-b +\sqrt {-4 a c +b^{2}}\right ) c}}\right )}{a^{2}}\)	\(244\)
risch	\(\text {Expression too large to display}\)	\(1298\)

[In]

int((f*x^4+e*x^2+d)/x^4/(c*x^4+b*x^2+a),x,method=_RETURNVERBOSE)

[Out]

-1/3*d/a/x^3-(a*e-b*d)/a^2/x+4/a^2*c*(1/8*(-a*e*(-4*a*c+b^2)^(1/2)+b*d*(-4*a*c+b^2)^(1/2)-2*f*a^2+a*b*e+2*a*c*
d-b^2*d)/(-4*a*c+b^2)^(1/2)*2^(1/2)/((b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctan(c*x*2^(1/2)/((b+(-4*a*c+b^2)^(1/2)
)*c)^(1/2))-1/8*(-a*e*(-4*a*c+b^2)^(1/2)+b*d*(-4*a*c+b^2)^(1/2)+2*f*a^2-a*b*e-2*a*c*d+b^2*d)/(-4*a*c+b^2)^(1/2
)*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)*arctanh(c*x*2^(1/2)/((-b+(-4*a*c+b^2)^(1/2))*c)^(1/2)))

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 9850 vs. \(2 (226) = 452\).

Time = 11.66 (sec) , antiderivative size = 9850, normalized size of antiderivative = 36.89 \[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=\text {Too large to display} \]

[In]

integrate((f*x^4+e*x^2+d)/x^4/(c*x^4+b*x^2+a),x, algorithm="fricas")

[Out]

Too large to include

Sympy [F(-1)]

Timed out. \[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=\text {Timed out} \]

[In]

integrate((f*x**4+e*x**2+d)/x**4/(c*x**4+b*x**2+a),x)

[Out]

Timed out

Maxima [F]

\[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=\int { \frac {f x^{4} + e x^{2} + d}{{\left (c x^{4} + b x^{2} + a\right )} x^{4}} \,d x } \]

[In]

integrate((f*x^4+e*x^2+d)/x^4/(c*x^4+b*x^2+a),x, algorithm="maxima")

[Out]

-integrate((a*b*e - a^2*f - (b*c*d - a*c*e)*x^2 - (b^2 - a*c)*d)/(c*x^4 + b*x^2 + a), x)/a^2 + 1/3*(3*(b*d - a
*e)*x^2 - a*d)/(a^2*x^3)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 3804 vs. \(2 (226) = 452\).

Time = 0.99 (sec) , antiderivative size = 3804, normalized size of antiderivative = 14.25 \[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=\text {Too large to display} \]

[In]

integrate((f*x^4+e*x^2+d)/x^4/(c*x^4+b*x^2+a),x, algorithm="giac")

[Out]

1/4*((sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*b^6 - 9*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^4*c - 2*sqrt
(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*b^5*c - 2*b^6*c + 24*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c^2 +
10*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^3*c^2 + sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*b^4*c^2 + 18*a
*b^4*c^2 - 2*b^5*c^2 - 16*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^3*c^3 - 8*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*
a*c)*c)*a^2*b*c^3 - 5*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^2*c^3 - 48*a^2*b^2*c^3 + 14*a*b^3*c^3 + 4*sq
rt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*c^4 + 32*a^3*c^4 - 24*a^2*b*c^4 + sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c
+ sqrt(b^2 - 4*a*c)*c)*b^5 - 7*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^3*c - 2*sqrt(2)*
sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*b^4*c + 12*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4
*a*c)*c)*a^2*b*c^2 + 6*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^2*c^2 + sqrt(2)*sqrt(b^2
- 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*b^3*c^2 - 3*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)
*a*b*c^3 + 2*(b^2 - 4*a*c)*b^4*c - 10*(b^2 - 4*a*c)*a*b^2*c^2 + 2*(b^2 - 4*a*c)*b^3*c^2 + 8*(b^2 - 4*a*c)*a^2*
c^3 - 6*(b^2 - 4*a*c)*a*b*c^3)*d - (sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^5 - 8*sqrt(2)*sqrt(b*c + sqrt(
b^2 - 4*a*c)*c)*a^2*b^3*c - 2*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^4*c - 2*a*b^5*c + 16*sqrt(2)*sqrt(b*
c + sqrt(b^2 - 4*a*c)*c)*a^3*b*c^2 + 8*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c^2 + sqrt(2)*sqrt(b*c
+ sqrt(b^2 - 4*a*c)*c)*a*b^3*c^2 + 16*a^2*b^3*c^2 - 2*a*b^4*c^2 - 4*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^
2*b*c^3 - 32*a^3*b*c^3 + 12*a^2*b^2*c^3 - 16*a^3*c^4 + sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*
c)*a*b^4 - 6*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c - 2*sqrt(2)*sqrt(b^2 - 4*a*c)
*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a*b^3*c + 8*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^3*c^2
+ 4*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b*c^2 + sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c
+ sqrt(b^2 - 4*a*c)*c)*a*b^2*c^2 - 2*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*c^3 + 2*(b^
2 - 4*a*c)*a*b^3*c - 8*(b^2 - 4*a*c)*a^2*b*c^2 + 2*(b^2 - 4*a*c)*a*b^2*c^2 - 4*(b^2 - 4*a*c)*a^2*c^3)*e + (sqr
t(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^4 - 8*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^3*b^2*c - 2*sqrt(2)
*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^3*c - 2*a^2*b^4*c + 16*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^4*c^2
+ 8*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^3*b*c^2 + sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c^2 +
16*a^3*b^2*c^2 - 2*a^2*b^3*c^2 - 4*sqrt(2)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^3*c^3 - 32*a^4*c^3 + 8*a^3*b*c^3
+ sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^3 - 4*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + s
qrt(b^2 - 4*a*c)*c)*a^3*b*c - 2*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c + sqrt(2)*
sqrt(b^2 - 4*a*c)*sqrt(b*c + sqrt(b^2 - 4*a*c)*c)*a^2*b*c^2 + 2*(b^2 - 4*a*c)*a^2*b^2*c - 8*(b^2 - 4*a*c)*a^3*
c^2 + 2*(b^2 - 4*a*c)*a^2*b*c^2)*f)*arctan(2*sqrt(1/2)*x/sqrt((a^2*b + sqrt(a^4*b^2 - 4*a^5*c))/(a^2*c)))/((a^
3*b^4 - 8*a^4*b^2*c - 2*a^3*b^3*c + 16*a^5*c^2 + 8*a^4*b*c^2 + a^3*b^2*c^2 - 4*a^4*c^3)*abs(c)) + 1/4*((sqrt(2
)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*b^6 - 9*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^4*c - 2*sqrt(2)*sqrt(b*c
- sqrt(b^2 - 4*a*c)*c)*b^5*c + 2*b^6*c + 24*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c^2 + 10*sqrt(2)*
sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^3*c^2 + sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*b^4*c^2 - 18*a*b^4*c^2 - 2
*b^5*c^2 - 16*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^3*c^3 - 8*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*
b*c^3 - 5*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^2*c^3 + 48*a^2*b^2*c^3 + 14*a*b^3*c^3 + 4*sqrt(2)*sqrt(b
*c - sqrt(b^2 - 4*a*c)*c)*a^2*c^4 - 32*a^3*c^4 - 24*a^2*b*c^4 + sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2
- 4*a*c)*c)*b^5 - 7*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^3*c - 2*sqrt(2)*sqrt(b^2 - 4
*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*b^4*c + 12*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2
*b*c^2 + 6*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^2*c^2 + sqrt(2)*sqrt(b^2 - 4*a*c)*sqr
t(b*c - sqrt(b^2 - 4*a*c)*c)*b^3*c^2 - 3*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b*c^3 - 2
*(b^2 - 4*a*c)*b^4*c + 10*(b^2 - 4*a*c)*a*b^2*c^2 + 2*(b^2 - 4*a*c)*b^3*c^2 - 8*(b^2 - 4*a*c)*a^2*c^3 - 6*(b^2
- 4*a*c)*a*b*c^3)*d - (sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^5 - 8*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)
*c)*a^2*b^3*c - 2*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^4*c + 2*a*b^5*c + 16*sqrt(2)*sqrt(b*c - sqrt(b^2
- 4*a*c)*c)*a^3*b*c^2 + 8*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c^2 + sqrt(2)*sqrt(b*c - sqrt(b^2 -
4*a*c)*c)*a*b^3*c^2 - 16*a^2*b^3*c^2 - 2*a*b^4*c^2 - 4*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b*c^3 + 32
*a^3*b*c^3 + 12*a^2*b^2*c^3 - 16*a^3*c^4 + sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a*b^4 - 6
*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c - 2*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c -
sqrt(b^2 - 4*a*c)*c)*a*b^3*c + 8*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^3*c^2 + 4*sqrt(2)
*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b*c^2 + sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 -
4*a*c)*c)*a*b^2*c^2 - 2*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*c^3 - 2*(b^2 - 4*a*c)*a
*b^3*c + 8*(b^2 - 4*a*c)*a^2*b*c^2 + 2*(b^2 - 4*a*c)*a*b^2*c^2 - 4*(b^2 - 4*a*c)*a^2*c^3)*e + (sqrt(2)*sqrt(b*
c - sqrt(b^2 - 4*a*c)*c)*a^2*b^4 - 8*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^3*b^2*c - 2*sqrt(2)*sqrt(b*c -
sqrt(b^2 - 4*a*c)*c)*a^2*b^3*c + 2*a^2*b^4*c + 16*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^4*c^2 + 8*sqrt(2)*
sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^3*b*c^2 + sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c^2 - 16*a^3*b^2*c
^2 - 2*a^2*b^3*c^2 - 4*sqrt(2)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^3*c^3 + 32*a^4*c^3 + 8*a^3*b*c^3 + sqrt(2)*sq
rt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b^3 - 4*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*
a*c)*c)*a^3*b*c - 2*sqrt(2)*sqrt(b^2 - 4*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b^2*c + sqrt(2)*sqrt(b^2 - 4
*a*c)*sqrt(b*c - sqrt(b^2 - 4*a*c)*c)*a^2*b*c^2 - 2*(b^2 - 4*a*c)*a^2*b^2*c + 8*(b^2 - 4*a*c)*a^3*c^2 + 2*(b^2
- 4*a*c)*a^2*b*c^2)*f)*arctan(2*sqrt(1/2)*x/sqrt((a^2*b - sqrt(a^4*b^2 - 4*a^5*c))/(a^2*c)))/((a^3*b^4 - 8*a^
4*b^2*c - 2*a^3*b^3*c + 16*a^5*c^2 + 8*a^4*b*c^2 + a^3*b^2*c^2 - 4*a^4*c^3)*abs(c)) + 1/3*(3*b*d*x^2 - 3*a*e*x
^2 - a*d)/(a^2*x^3)

Mupad [B] (verification not implemented)

Time = 10.74 (sec) , antiderivative size = 15505, normalized size of antiderivative = 58.07 \[ \int \frac {d+e x^2+f x^4}{x^4 \left (a+b x^2+c x^4\right )} \, dx=\text {Too large to display} \]

[In]

int((d + e*x^2 + f*x^4)/(x^4*(a + b*x^2 + c*x^4)),x)

[Out]

atan(((x*(4*a^8*c^5*d^2 - 4*a^9*c^4*e^2 + 4*a^10*c^3*f^2 + 2*a^6*b^4*c^3*d^2 - 8*a^7*b^2*c^4*d^2 + 2*a^8*b^2*c
^3*e^2 - 8*a^9*c^4*d*f + 12*a^8*b*c^4*d*e - 4*a^9*b*c^3*e*f - 4*a^7*b^3*c^3*d*e + 4*a^8*b^2*c^3*d*f) - (-(b^7*
d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3
*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^
3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^
5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(
1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c
*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e
+ 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 -
8*a^6*b^2*c)))^(1/2)*(x*(32*a^11*b*c^3 - 8*a^10*b^3*c^2)*(-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3
)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e
^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1
/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e -
2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f +
24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*
e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4
*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) - 16*a^10*c^4*d + 16*a^
11*c^3*f - 4*a^8*b^4*c^2*d + 20*a^9*b^2*c^3*d + 4*a^9*b^3*c^2*e - 4*a^10*b^2*c^2*f - 16*a^10*b*c^3*e))*(-(b^7*
d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3
*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^
3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^
5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(
1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c
*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e
+ 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 -
8*a^6*b^2*c)))^(1/2)*1i + (x*(4*a^8*c^5*d^2 - 4*a^9*c^4*e^2 + 4*a^10*c^3*f^2 + 2*a^6*b^4*c^3*d^2 - 8*a^7*b^2*
c^4*d^2 + 2*a^8*b^2*c^3*e^2 - 8*a^9*c^4*d*f + 12*a^8*b*c^4*d*e - 4*a^9*b*c^3*e*f - 4*a^7*b^3*c^3*d*e + 4*a^8*b
^2*c^3*d*f) - (-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b
^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*
a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) -
9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e
*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^
2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2)
- 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*(x*(32*a^11*b*c^3 - 8*a^10*b^3*c^2)*(-(b^7*d^2 + a^2*b^5*e^2 + b^4*
d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c
*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*
(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*
f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e
- 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^
(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*
c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) +
16*a^10*c^4*d - 16*a^11*c^3*f + 4*a^8*b^4*c^2*d - 20*a^9*b^2*c^3*d - 4*a^9*b^3*c^2*e + 4*a^10*b^2*c^2*f + 16*a
^10*b*c^3*e))*(-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b
^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*
a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) -
9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e
*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^
2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2)
- 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a
^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*1i)/((x*(4*a^8*c^5*d^2 - 4*a^9*c^4*e^2 + 4*a^10*c^3*f^2 + 2*a^6*b^4
*c^3*d^2 - 8*a^7*b^2*c^4*d^2 + 2*a^8*b^2*c^3*e^2 - 8*a^9*c^4*d*f + 12*a^8*b*c^4*d*e - 4*a^9*b*c^3*e*f - 4*a^7*
b^3*c^3*d*e + 4*a^8*b^2*c^3*d*f) - (-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 +
a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c
- b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*
a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*
c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^
3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(
4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c -
b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*(x*(32*a^11*b*c^3 - 8*a^10*b^3*c^2)*(-(b^7*d^2
+ a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*
c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c
^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b
*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2
) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*
f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e +
2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*
a^6*b^2*c)))^(1/2) - 16*a^10*c^4*d + 16*a^11*c^3*f - 4*a^8*b^4*c^2*d + 20*a^9*b^2*c^3*d + 4*a^9*b^3*c^2*e - 4*
a^10*b^2*c^2*f - 16*a^10*b*c^3*e))*(-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 +
a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c
- b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*
a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*
c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^
3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(
4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c -
b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) - (x*(4*a^8*c^5*d^2 - 4*a^9*c^4*e^2 + 4*a^10*c
^3*f^2 + 2*a^6*b^4*c^3*d^2 - 8*a^7*b^2*c^4*d^2 + 2*a^8*b^2*c^3*e^2 - 8*a^9*c^4*d*f + 12*a^8*b*c^4*d*e - 4*a^9*
b*c^3*e*f - 4*a^7*b^3*c^3*d*e + 4*a^8*b^2*c^3*d*f) - (-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/
2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 -
a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) +
a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3
*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^
4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f -
3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*
b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*(x*(32*a^11*b*c^3 - 8*a^10*b
^3*c^2)*(-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)
^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*
d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b
^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*
a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^
(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*
a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4
+ 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) + 16*a^10*c^4*d - 16*a^11*c^3*f + 4*a^8*b^4*c^2*d - 20*a^9*b^2*c^3*d - 4*
a^9*b^3*c^2*e + 4*a^10*b^2*c^2*f + 16*a^10*b*c^3*e))*(-(b^7*d^2 + a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/
2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 -
a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) +
a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3
*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^
4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f -
3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e + 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*
b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) + 2*a^8*c^4*e^3 - 2*a^6*b*c^
5*d^3 + 2*a^7*c^5*d^2*e + 2*a^9*c^3*e*f^2 - 4*a^8*c^4*d*e*f - 4*a^7*b*c^4*d*e^2 + 4*a^7*b*c^4*d^2*f - 2*a^8*b*
c^3*d*f^2 - 2*a^8*b*c^3*e^2*f + 2*a^6*b^2*c^4*d^2*e - 2*a^6*b^3*c^3*d^2*f + 4*a^7*b^2*c^3*d*e*f))*(-(b^7*d^2 +
a^2*b^5*e^2 + b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 + a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^
3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 - a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2
*d^2 + a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) + a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c
*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f - 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2)
+ 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f - 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) - 2*a^3*c*d*f*
(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f - 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e + 2*
a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) + 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^
6*b^2*c)))^(1/2)*2i - (d/(3*a) + (x^2*(a*e - b*d))/a^2)/x^3 + atan(((x*(4*a^8*c^5*d^2 - 4*a^9*c^4*e^2 + 4*a^10
*c^3*f^2 + 2*a^6*b^4*c^3*d^2 - 8*a^7*b^2*c^4*d^2 + 2*a^8*b^2*c^3*e^2 - 8*a^9*c^4*d*f + 12*a^8*b*c^4*d*e - 4*a^
9*b*c^3*e*f - 4*a^7*b^3*c^3*d*e + 4*a^8*b^2*c^3*d*f) - (-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(
1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2
+ a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2)
- a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a
^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*
a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f
+ 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^
2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*(x*(32*a^11*b*c^3 - 8*a^10
*b^3*c^2)*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^
3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^
6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a
*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(
4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3
)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 3
6*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b
^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) - 16*a^10*c^4*d + 16*a^11*c^3*f - 4*a^8*b^4*c^2*d + 20*a^9*b^2*c^3*d +
4*a^9*b^3*c^2*e - 4*a^10*b^2*c^2*f - 16*a^10*b*c^3*e))*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(
1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2
+ a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2)
- a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a
^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*
a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f
+ 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^
2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*1i + (x*(4*a^8*c^5*d^2 - 4
*a^9*c^4*e^2 + 4*a^10*c^3*f^2 + 2*a^6*b^4*c^3*d^2 - 8*a^7*b^2*c^4*d^2 + 2*a^8*b^2*c^3*e^2 - 8*a^9*c^4*d*f + 12
*a^8*b*c^4*d*e - 4*a^9*b*c^3*e*f - 4*a^7*b^3*c^3*d*e + 4*a^8*b^2*c^3*d*f) - (-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2
*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^
2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(
4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f +
16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 1
4*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/
2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c -
b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*(x*(3
2*a^11*b*c^3 - 8*a^10*b^3*c^2)*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4
*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b
^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c
- b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*
e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*
e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*
c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2
)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) + 16*a^10*c^4*d - 16*a^11*c^3*f + 4*a^8*b^4*c^2*d
- 20*a^9*b^2*c^3*d - 4*a^9*b^3*c^2*e + 4*a^10*b^2*c^2*f + 16*a^10*b*c^3*e))*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2
*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^
2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(
4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f +
16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 1
4*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/
2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c -
b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*1i)/(
(x*(4*a^8*c^5*d^2 - 4*a^9*c^4*e^2 + 4*a^10*c^3*f^2 + 2*a^6*b^4*c^3*d^2 - 8*a^7*b^2*c^4*d^2 + 2*a^8*b^2*c^3*e^2
- 8*a^9*c^4*d*f + 12*a^8*b*c^4*d*e - 4*a^9*b*c^3*e*f - 4*a^7*b^3*c^3*d*e + 4*a^8*b^2*c^3*d*f) - (-(b^7*d^2 +
a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3
*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*
d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*
f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) +
16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(
-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a
^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6
*b^2*c)))^(1/2)*(x*(32*a^11*b*c^3 - 8*a^10*b^3*c^2)*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2
) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a
^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) -
a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*
b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4
*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f +
3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b
*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) - 16*a^10*c^4*d + 16*a^11*c^3
*f - 4*a^8*b^4*c^2*d + 20*a^9*b^2*c^3*d + 4*a^9*b^3*c^2*e - 4*a^10*b^2*c^2*f - 16*a^10*b*c^3*e))*(-(b^7*d^2 +
a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3
*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*
d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*
f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) +
16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(
-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a
^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6
*b^2*c)))^(1/2) - (x*(4*a^8*c^5*d^2 - 4*a^9*c^4*e^2 + 4*a^10*c^3*f^2 + 2*a^6*b^4*c^3*d^2 - 8*a^7*b^2*c^4*d^2 +
2*a^8*b^2*c^3*e^2 - 8*a^9*c^4*d*f + 12*a^8*b*c^4*d*e - 4*a^9*b*c^3*e*f - 4*a^7*b^3*c^3*d*e + 4*a^8*b^2*c^3*d*
f) - (-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1
/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e
+ 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*
c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c
- b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/
2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3
*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 +
16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*(x*(32*a^11*b*c^3 - 8*a^10*b^3*c^2)*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*
a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12
*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c
- b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^
4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*
b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 1
2*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^
3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) + 16*a^10*c
^4*d - 16*a^11*c^3*f + 4*a^8*b^4*c^2*d - 20*a^9*b^2*c^3*d - 4*a^9*b^3*c^2*e + 4*a^10*b^2*c^2*f + 16*a^10*b*c^3
*e))*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2) + a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1
/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a^3*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e
+ 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) - a^2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*
c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c
- b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/
2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f + 3*a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3
*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b*c*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 +
16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2) + 2*a^8*c^4*e^3 - 2*a^6*b*c^5*d^3 + 2*a^7*c^5*d^2*e + 2*a^9*c^3*e*f^2 - 4*a^
8*c^4*d*e*f - 4*a^7*b*c^4*d*e^2 + 4*a^7*b*c^4*d^2*f - 2*a^8*b*c^3*d*f^2 - 2*a^8*b*c^3*e^2*f + 2*a^6*b^2*c^4*d^
2*e - 2*a^6*b^3*c^3*d^2*f + 4*a^7*b^2*c^3*d*e*f))*(-(b^7*d^2 + a^2*b^5*e^2 - b^4*d^2*(-(4*a*c - b^2)^3)^(1/2)
+ a^4*b^3*f^2 - a^4*f^2*(-(4*a*c - b^2)^3)^(1/2) - 20*a^3*b*c^3*d^2 - 7*a^3*b^3*c*e^2 + 12*a^4*b*c^2*e^2 + a^3
*c*e^2*(-(4*a*c - b^2)^3)^(1/2) - 2*a*b^6*d*e + 25*a^2*b^3*c^2*d^2 - a^2*b^2*e^2*(-(4*a*c - b^2)^3)^(1/2) - a^
2*c^2*d^2*(-(4*a*c - b^2)^3)^(1/2) - 9*a*b^5*c*d^2 - 4*a^5*b*c*f^2 + 2*a^2*b^5*d*f + 16*a^4*c^3*d*e - 2*a^3*b^
4*e*f - 16*a^5*c^2*e*f + 2*a*b^3*d*e*(-(4*a*c - b^2)^3)^(1/2) + 16*a^2*b^4*c*d*e - 14*a^3*b^3*c*d*f + 24*a^4*b
*c^2*d*f + 2*a^3*b*e*f*(-(4*a*c - b^2)^3)^(1/2) + 2*a^3*c*d*f*(-(4*a*c - b^2)^3)^(1/2) + 12*a^4*b^2*c*e*f + 3*
a*b^2*c*d^2*(-(4*a*c - b^2)^3)^(1/2) - 36*a^3*b^2*c^2*d*e - 2*a^2*b^2*d*f*(-(4*a*c - b^2)^3)^(1/2) - 4*a^2*b*c
*d*e*(-(4*a*c - b^2)^3)^(1/2))/(8*(a^5*b^4 + 16*a^7*c^2 - 8*a^6*b^2*c)))^(1/2)*2i